Optimal reaction time for surface-mediated diffusion 
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We present an exact calculation of the mean first-passage time to a small target on the surface of 
a 2D or 3D spherical domain, for a molecule performing surface-mediated diffusion. This minimal 
model of interfacial reactions, which explicitly takes into account the combination of surface and bulk 
diffusion, shows the importance of correlations induced by the coupling of the switching dynamics 
to the geometry of the confinement, ignored so far. Interestingly, our results show that, in the 
context of interfacial systems in confinement, the reaction time can be minimized as a function of 
the desorption rate from the surface, which puts forward a general mechanism of enhancement and 
regulation of chemical and biological reactivity. 

PACS numbers: 



Among reactions limited by transport, interfacial reac- 
tions, for which molecules react on target sites located on 
the surface of the confining domain [1] , play an important 
role in situations as various as heterogeneous catalysis 
[2J, reactions in porous media or biochemical reactions 
on DNA [3] and in vesicular systems fH-E] . Besides being 
a problem of great practical importance, the modeling of 
such systems raises two types of theoretical issues: (i) to 
determine the impact of geometrical parameters of the 
confining domain, such as the volume, on reaction kinet- 
ics, and (ii) to account for the reactive trajectories which 
combine bulk and surface-mediated diffusion phases due 
to the affinity of the molecules for the surface (see fig.l). 

The point (i) has been studied in the case of a perfectly 
reflecting surface [SHTU] and a universal scaling with the 
confining volume was found for the mean first-passage 
time (MFPT) to a reactive target PUCE]- On the other 
hand, reactive paths of point (ii) can be described as tra- 
jectories involving two states, one adsorbed on the surface 
and one desorbed in the bulk. Such two-state paths have 
been studied in the broader context of intermittent search 
strategies under the hypothesis that the times spent in 
each state are controlled by an internal clock independent 
of any geometrical parameter [T3TfT5] . In most cases, the 
sojourn times in each state have been assumed to be ex- 
ponentially distributed Q3], with the notable exception 
of Levy [Hj and deterministic laws [TTJ, [TH] . 

However, in the case of interfacial reactions in confine- 
ment, the time spent in a bulk excursion is controlled 
by the statistics of return to the surface and therefore 
by the geometry of the confining domain [TWI22J . Hence, 
this return time is not an external parameter but is gen- 
erated by the very dynamics of the diffusing molecule in 
confinement. As a result, the usual methods to calcu- 
late mean search times for intermittent processes only 
provide a mean-field (MF) approximation of the reaction 
time which completely ignores the correlations induced 
by the coupling of the switching dynamics to the geom- 
etry of the confinement (see |14j for review, and more 



recently [2"3"]). 

Here we calculate exactly in the representative exam- 
ple of a spherical confining domain the MFPT to a re- 
active site on the sphere for a Brownian molecule alter- 
nating phases of bulk and surface diffusion. Our ana- 
lytical approach fully complies with points (i) and (ii) 
and shows that correlations actually strongly impact on 
reaction times, which are substantially underestimated 
by standard MF treatments. In addition, we discuss the 
problem of the minimization of the reaction time with re- 
spect to the mean adsorption time of the molecule on the 
surface, which is a priori not clear. Indeed, the benefit of 
bulk diffusion, even if much faster than surface diffusion, 
is questionable, since the mean time spent in bulk excur- 
sions diverges with the volume of the confining domain. 
Surprisingly enough, we will show that, even for bulk and 
surface diffusion coefficients of the same order of magni- 
tude, the reaction time can be minimized, whereas the 
MF treatment of [23] predicts a monotonic behavior. 

At the theoretical level, these results bring an exact 
solution to an extension to the so-called narrow escape 
problem (time needed to escape through a small window 
of an otherwise reflecting domain) , which has attracted a 
lot of attention in last few years both in the mathematical 
El El ES] and physical (23 ES] literature, partly due 
to the challenge of taking into account mixed boundary 
conditions. At the level of bio and chemical physics, it 
puts forward a general mechanism of enhancement and 
regulation of chemical reactivity by tuning the affinity of 
the rcactants for the surface of the confining domain. 

As an archetype of confined interfacial systems, we 
consider a molecule in a spherical domain S of radius R 
(see fig. [I]), alternating phases of surface diffusion on dS 
with diffusion coefficient D\ and phases of bulk diffusion 
in S with diffusion coefficient D 2 ■ The time spent during 
each surface phase is assumed to follow an exponential 
law with desorption rate A, which is reminiscent of a first- 
order kinetics. At each desorption event, the molecule 
is assumed to be radially ejected at a distance a from 



the surface (otherwise it is instantaneously readsorbed) . 
Although formulated for any value of this parameter a 
smaller than R, in most situations of physical interest 
a -c R. For the sake of simplicity, we first present the 
case of a point-like target on a 2D sphere (disk), and then 
generalize our approach to targets of angular extension 
2e on 2D and 3D spheres. 

target 




FIG. 1: Model of surface-mediated reaction in confinement. 

Point-like target in 2D. We start with the example of 
a point-like target (e = 0) on the surface of a 2D sphere. 
Note that, as opposed to narrow-escape problems for 
one-state diffusion [6], the limit e — > is not singular 
and the MFPT at the target is finite, thanks to the one- 
dimensional surface diffusion. This MFPT satisfies the 
following backward equations [27]: 



D 2 A s t 2 (r, 



-1, 



(2) 



where t\ stands for the MFPT starting from the ad- 
sorbed state at a position defined on the surface by 
the polar angle 9, and t 2 for the MFPT starting from 
the point (r, 9) in the bulk. Here, Ags = d 2 /R 2 and 
As = d 2 + d r /r + d 2 /r 2 . In Eqs@,([2) the first term of 
the lhs accounts for the diffusion respectively on the sur- 
face and in the bulk, while the second term of Eq. de- 
scribes desorption events. They have to be completed by 
boundary conditions: t 2 (R,9) = ti(9), which describes 
the adsorption events and t\{9 = 0) = = t\{6 = 2tt), 
which expresses that the target is an absorbing point in 
the problem. As we proceed to show, these equations can 
be solved exactly. 

Considering t 2 as a source term in the Poisson type 
equation (fTl) with absorbing conditions at 9 = and 9 — 
2ir, whose Green function is well known [28 , t\ writes: 
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whe re we have used the dimensionless variable u = 
-Ra/A/Di and the notations 9 < = min(0, 9') and 9 > — 
max((9, 9'). On the other hand, Eq.pJ) is easily shown to 
be satisfied by the following Fourier series: 
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where the unknown coefficients a n have to be calcu- 
lated. We aim at determining the reaction time, de- 
fined here as the MFPT at the target, for a molecule 
initially uniformly distributed on the circumference, ie 
(ti) = ± J^t 1 (9)d9. Substituting Eqs (J3|) , Q in the 
boundary conditions leads after straightforwara integra- 
tions to 



to 
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where x = 1 — a/R. The projection of Eq([5]) on the func- 
tions cos(n#) leads to an infinite hierarchy of equations 
for a n , n £ N, which can be decoupled. It finally yields 
an exact expression of the reaction time: 
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Several comments are in order, (i) As expected, both 
limits A — > and a — > of the reaction time are given 
by tt 2 R 2 /(3Di) which corresponds to a pure ID diffusion, 
with a molecule initially uniformly distributed in [0, 2-7ri?] 
with absorbing boundaries [37]. (ii) The limit a —*■ R of 
Eq.(|6]) corresponds actually to rebinding positions uni- 
formly distributed on the surface, and can therefore be 
recovered from the MF treatment of [14], which ignores 
the correlations between starting and ending points of 
bulk excursions. Actually, it can be checked that this 
type of MF treatment, which has been successfully used 
in the context of target search on DNA [14] , substantially 
underestimates the reaction time in the physical regime 
a <C R in the present case of interfacial reactions, where 
correlations play a crucial role. Note that the MF ap- 
proach of|23j even misses the non-monotonicity with A. 
(iii) Eq.([6]) has a clear physical interpretation: the re- 
action time is the product of the mean duration of each 
elementary step composed of one surface exploration and 
one excursion in the bulk phase, by the mean number of 
excursions before reaction. This factorized structure is 
expected in the limit of large number of excursions since 



the durations of elementary steps are independent vari- 
ables, all independent of the number of excursions, except 
for the very last one. 

Importantly, a first-order small A expansion of the re- 
action time allows one to show that the reaction time is 
a non-monotonic function of the desorption rate A if 
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which becomes D 2 jD x > vr 2 /(12C(3)) « 0.68... in the 
physically relevant limit a/R <C 1 (ie x — > 1), where £ 
is the Riemann function. In other words, even if D 2 is 
smaller than D\ , bulk excursions turn out to speed up the 
reaction. Focusing on this limit a/R <C 1 and D\ <C D 2 , 
a detailed analysis of the reaction time given explicitly 
by Eq.dFj]) shows that the optimal A is given to leading 
order by A opt « 2D 2 ln(2D 2 /D 1 )/(aR). In turn, the gain, 
defined as the ratio of the reaction time in absence of bulk 
excursion over the optimal reaction time, is found to be 
proportional to D 2 /Di, up to logarithmic corrections. 
This shows that for fast bulk diffusion, reactivity can be 
significantly enhanced by surface-mediated diffusion by 
a proper tuning of the affinity of the reactants for the 
surface (see fig[2] for e = 0) . 
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FIG. 2: Mean reaction time as a function of the desorption 
rate in the 2D case, with D\ — 1, D2 = 5, a = 0.1 and R = 1 
(in arbitrary units): analytical perturbative expression Eq.® 
vs Monte-Carlo numerical simulations for different target sizes 
e. Insert: mean reaction time as a function of the target size 
for a fixed desorption rate A = 0.125. 



Extended target in 2D. The previous analysis can be 
generalized to the important case of an extended tar- 
get zone (of angular extension 2e), especially relevant to 
the case of escape problems [6 . The calculation leads 
in this case to an infinite hierarchy of coupled equations 
for a n , n g N, in contrast to the point-like target case. 
This additional complexity results from the fact that the 
target can now be reached not only from the surface, but 
also directly from the bulk. The first terms of an exact 



perturbative expansion in e can be obtained and read: 
(t 1 (e)) = <t 1 (e = 0)}+ W 2 



A + i^ (1 -^ 



i+ 5:^7i 



2w 2 (l - x m ) 
(1 - x m ) + m 2 
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Note that the coefficients of e k of this expansion diverge 
with uj, so that in practice, the smaller u>, the wider the 
range of applicability in e. Fig[2]shows an excellent quan- 
titative agreement even for rather large values of e. 

In addition, the benefit of bulk excursions can still be 
analyzed with e 7^ 0. A small A expansion of the reaction 
time can be worked out, and shows that bulk excursions 
reduce the reaction time provided that: 
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FIG. 3: Mean reaction time as a function of the desorption 
rate in the 3D case with Di — 1, D2 — 50, a = 0.02 and 
R — 1 (in arbitrary units): analytical perturbative expression 
Eq.( |10[ | vs Monte-Carlo simulations for a target of size e = 
0.02. Insert: mean reaction time as a function of the target 
size for a fixed desorption rate A = 0.125. 

3D case. The previous analysis can also be adapted 
to the case of a confining 3D sphere, relevant to many 
situations such as reactions in micellar or vesicular sys- 
tems [3 |B]. EqsQ-^ are still valid in this case with 
A-ds — de(sin9de)/(R 2 sin 8) and A5 is the usual 3D 
Laplacian. Along the same lines (technical details are 
given in Supplementary Information), the first terms of 
an exact small e expansion can be obtained and read : 
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Similarly to the 2D case, one can show that bulk excur- 



sions are beneficial under the condition that now reads 
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where6„ = ( (n cos e + n + 1 )P„ (cos e) + P n _i (cose)) 2 and 
P n stand for Legendre polynomials, and that the reaction 
time can be minimized. Fig|3] shows a good quantitative 
agreement of the expansion Eq.(lO) with Monte-Carlo 



simulations and confirms that the reaction time can be 
decreased by bulk excursions, as in the 2D situation. 



Entire FPT distribution. Last, using the results given 
recently in :29, [30] , we stress that an estimate of the en- 
tire FPT distribution of the reduced variable ii/(ti) can 
be inferred in the large size limit R — > oo from the knowl- 
edge of the MFPT calculated above. While an exponen- 
tial form is predicted in the 3D case (since the exploration 
of the boundary is marginally compact), the FPT distri- 
bution in the 2D case involves several time scales, and is 
given by a sum of exponentials (see Eq.(5) of [30], case 
of compact exploration) . Fig. H shows a good agreement 
of this theoretical prediction with numerical simulations 
in both cases. 




FIG. 4: FPT distribution of the reduced variable ti/(ti) : 
theory (plain line) vs simulations (symbols). In the 2D case, 
Dx = 1, D 2 = 10, a = 0.1, e = 0.01 and R = 1 (in arbitrary 
units). Insert: 3D case, with D\ — 1, D2 = 50, a = 0.02, 
A = 70, e = 0.02 and R = 1. 



In conclusion, we have presented an exact calcula- 
tion of the mean first-passage time to a small target 
on the surface of a 2D and 3D spherical domain, for 
a molecule performing surface-mediated diffusion. This 
minimal model, which explicitly takes into account the 
combination of surface and bulk diffusion, shows that 
bulk excursions can speed up the reaction as they re- 
duce oversampling of the boundary and underlines the 
importance of correlations induced by the coupling of the 
switching dynamics to the geometry of the confinement. 
In the context of intcrfacial systems in confinement, our 
results show that the reaction time can be minimized as 
a function of the desorption rate from the surface, which 
puts forward a general mechanism of enhancement and 
regulation of chemical reactivity. 
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